arXiv:1501.01078v3 [cond-mat.mtrl-sci] 31 Mar 2015 


Dispersion relation in chiral media : 
Credibility of Drude-Born-Fedorov equations 

Kikuo Cho 

Institute of Laser Engineering, Osaka University, Suita 565-0871, Osaka Japan 


Abstraet 

Disperion relation of EM field in a ehiral medium is diseussed from the viewpoint of 
eonstitutive equations to be used as a partner of Maxwell equations. The popular form 
of Drude-Born-Eedorov (DBE) eonstitutive equations is eritieized via a eomparison with 
the first-prineiples maeroseopie eonstitutive equations. The two sets of equations show 
a deeisive differenee in the dispersion eurve in the resonant region of ehiral, left-handed 
eharaeter, in the form of presenee or absenee of linear erossing at k=0. DBE equations 
eould be used at most only as a phenomenology in off-resonant region, while the first- 
prineiples ones ean be used for both phenomenologieal and mieroseopie analyses. 


1 Introduction 

Symmetry plays an important role in the eleetromagnetie (EM) response of matter. It 
is revealed in the form of suseeptibilities relating eleetrie and magnetie polarization 
(P and M) with souree EM field. In high symmetry ease, P and M eonsist of (the 
superpositions of) independent groups of exeitations belonging to different irredueible 
representations of the symmetry group in eonsideration. This allows us to treat eleetrie 
and magnetie properties of matter independently. When a medium laeks in eertain mir¬ 
ror symmetry, i.e., the ease of ehiral symmetry, however, some (or all the) eomponents 
of P and M eannot be distinguished, so that they ean be indueed by both eleetrie and 
magnetie souree fields. In addition, there is also a mixing between eleetrie dipole (El) 
and eleetrie quadrupole (E2) transitions. 

The study of ehiral symmetry in the EM response of matter has a long history (Intro- 
duetion of [[U). Chiral substanees have been eonsidered as uneonventional materials for 
a long time, but now it is regarded as an important souree of new materials and states, 
providing hot topies in the studies of metamaterials [|2l, multiferroies [l3]|, and supereon- 
duetivity [j4]|. 

In spite of its long history, theoretieal deseription of ehirality does not seem to be 
standardized. In the doeuments of lUPAP and lUPAC dealing with the standard defini¬ 
tions of physieal and ehemieal quantities [[510, there is no mentioning about the ehiral 
suseeptibilities. Correspondingly, there are two or more different forms of phenomeno- 
logieal eonstitutive equations in use for maeroseopie response. Though the effeet of 
ehiral symmetry is expeeted also in mieroseopie responses, its first-prineiples thoery has 
been made only very reeently 0. Prom the viewpoit that all the different forms of EM 
response theories should belong to a single hierarehy with logieal ranking, one should 
be able to ehoose the most appropriate form of the eonstitutive equations for the maero¬ 
seopie ehiral response on the basis of the mieroseopie theory. 

A typieal effeet of ehirality is the differenee in the phase veloeity of EM waves with 
right- and left-eireular polarizations, whieh appears in the off-resonant region of sus- 
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ceptibilities. However, this is not the only aspect of our interest in discussing chirality. 
In fact, the dispersion curves in the resonant region of susceptibility show a remarkable 
behavior, by which we can select the correct constitutive equations. 

Macroscopic EM response of matter is usually calculated by the combination of 
Maxwell and constitutive equations. The standard form of the latter is 

D = eE , B = nH (1) 

with the dielectric constant (permittivity) e and permeability fx. However, if the medium 
in consideration has chiral symmetry, these constitutive equations need to be generalized. 
A popular form of such an extention is 

D = e{E + /3VxE), (2) 

B = ^i{H + /3VxH), (3) 

which is called Drude-Born-Fedorov equations (DBF eqs) [[H13 . The parameter (chiral 
admittance) describes the chirality of the medium. This is a phnomenology for uniform 
and isotropic media. 

However, this is not the only way of generalization. From the viewpoint that the 
fundamental variables of EM field are E and B, both electric and magnetic polariza¬ 
tions P and M should consist both of the {E and .B}-induced components, so that the 
definition D = E + dvrP, H = B — AttM leads to the extension 

D = eE + i^B ^ (4) 

H = {1/fi)B + iT]E , (5) 

where the terms with ^ and rj take care of the chirality. For later convenience, let us call 
them chiral constitutive equations (ChC eqs). Though they are a result of phenomeno¬ 
logical consideration on the one hand, a first-principles calculation of macroscopic con¬ 
stitutive equations can be put also in this form on the other hand JTl . 

As to the difference or similarity of DBF and ChC eqs, there is a controversy. There 
have been arguments in the metamaterials community that DBF and ChC eqs are es¬ 
sentially same ifTOll . and also it is argued that the former can be derived from the latter 
by assuming the uniformity and isotropy of matter ||Tl(Sec.4.4). But there is also other 
group of people preferring ChC to DBF eqs. The purpose of this article is to show that 
there is a clear difference between the two, and that ChC eqs should be preferred. In 
view of the fact that DBF eqs are frequently used in metamaterials studies and also in 
recent textbook of standard electromagnetism |I3, it will be important to clarify the dif¬ 
ference between DBF and ChC eqs. 

We first note the relation between the parameters of DBF and ChC eqs. By means of 
the relation 

V X B = {iu/c)B , V X iT = {—iu/c)D , (6) 

DBF eqs can be rewritten as 


= eE + {iu/c)e(3B ^ 

= (ia;/c)e/9B-f (l//i)[l — (a;/9/c)^e/i]B . 
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If DBF and ChC eqs are equivalent, the DBF parameters can be written in terms of the 
ChC parameters by comparing (|71) and ([8]) with (jH) and ([5]) as 

e = e, ^ = T] = {u/c)€(3 , (l//i) = (l//i) - (a;/3/c)2e . (9) 

This relation will be shown later to lead to contradiction, which disproves the equiva¬ 
lence of DBF and ChC eqs. 


The first-principles derivation of micro- and macroscopic constitutive equations is 
done in the following way [|7]|. We assume a general form of non-relativistic Hamil¬ 
tonian (including relativistic correction terms, such as spin-orbit interaction and spin 
Zeeman term, etc.) for a many particle system in an EM field, and calculate the micro¬ 
scopic current density induced by the EM field, which is in general given as a functional 
of the transverse (T) part of vector potential and the longitudinal (E) external elec¬ 
tric field eI^I. The integral kernel of the functional is the microscopic susceptibility 
of a separable form with respect to position coordinates. When the relevant quantum 
mechanical states have spatial extension much less than the wavelength of the EM field, 
we may apply long wavelength approximation to the microscopic current density, which 
leads to the macroscopic constitutive equations to be used for macroscopic Maxwell eqs. 
In this macroscopic scheme, we need only a single 3x3 tensor to relate induced current 
density and source EM field, covering all the electric, magnetic and chiral polarizations 
of matter. This macroscopic constitutive equation is given in the form [|71 

J{k,uj) = Xem{k,uj) - {ic/uj)E[^l{k,uj)} . (10) 

The internal E field does not appear in the source field, since it is taken into account 
as the Coulomb potential in the matter Hamiltonian. The susceptibility Xem is written 
in terms of the quantum mechanical transition energies and the lower moments of the 
corresponding transition matrix elements of current density operator. 

Using the identity J = —iuP+ick x M in Eourier representation, we can rigorously 
rewrite the constitutive equation into the form 

P = XeEp + Xe^B , M = XhieE -f XhibB (11) 


The four susceptibilities XeE, XeB, XmE, XmB are again written in terms of the quantum 
mechanical transition energies and lower transition moments of electric dipole (El), 
electric quadrupole (E2), and magnetic dipole (Ml) characters. Details are given in 
sec.3.1 of im. The lowest order terms of them are 
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where V is the volume of a eell for periodie boundary eondition to define k, and and 
Mou are, respeetively, the El and Ml transition moments of eurrent density and (orbital 
and spin) magnetization operators between the matter eigenstates |0) (ground state) and 
\u) with transition energy E^q between them. (E2 moments appear in the terms 
in the next higher order.) Chiral symmetry allows the existenee of the transitions with 
mixed (El and Ml) or (El and E2) eharaeter, leading to the 0{k^) terms in Xem- 

Though there appear four suseeptibilities, the single suseeptibility nature is intaet, 
sinee the rewriting of (fTOl) into (fTTl) is reversible. Combining the new form of eonstitutive 
equations with the definition of D and H, we obtain 

D = E + AnP = {1 +A'KXeE)E + AnXeBB , (14) 

H = B - AttM = (1 - 47rXmB)-B - , (15) 

whieh is essentially equivalent to ChC eqs. The parameters of ChC eqs are given as 

e = 1 + dvrXeE , = dvr^es , iv = -4vrXmE , A = ^y -• (16) 

1 - dvrXmB 

all of whieh are tensors, with no assumption of isotropy and homogeneity as for DBE 
eqs. It should also be noted that the poles of XmB, i-e-, the magnetie transition energies, 
are, not the poles, but the zeros of fi. This is due to the definition of XmB, M = XmB-B as 
required in the first-prineiples approaeh, in eontrast to the eonventional one M = XmH- 
At this stage, the ChC eqs are not a phenomenology, but a first-prineiples theory. In 
eontrast, this kind of first-prineiples derivation does not exist for DBE eqs. 

2 Dispersion equation 

In order to show the differenee between DBE and ChC eqs, it is suffieient to give a single 
example. Eor this purpose,we eompare the dispersion relations obtained from DBE and 
ChC eqs. 

2.1 Case of DBF eqs 

If we solve DBE eqs and eq.® for V x iT and V x we obtain 

V X ET = aH + bE , 

V X E = dH + eE , 

where 

a = e = —e/i/9/A , b = —ice/uA , d = +icfx/uA , (19) 

and A = e/i/5^ — (?juP'. Erom eqs ®, ©, ® and V ■ .B = 0, both E and H are 
transverse, so that 

V X V X B = PE, V X V X Ef = pH (20) 

for Eourier eomponents. Then, by operating Vx to (fTTl) and (fTSl) . we obtain a set of 
homogeneous linear equations of EE, E. The eondition for the existenee of non-trivial 
solution gives us the dispersion equation 

detlfc^l-^2| ^ 0 , (21) 


(17) 

(18) 
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where ^ is a 2 x 2 matrix with the components a, b, d, e 


1-1 

b ' 

e 

(22) 

This dispersion equation can be rewritten as 



det A;1 -f ^ = 0 or 

det /cl — ^ = 0 , 

(23) 

so that the solution is 

k = ia i 

Vm , 

(24) 


with all the combinations of ± being allowed, which finally leads to a compact expres¬ 
sion 

— = ±-^ - C25) 

cj l±(u;/3/c)^' ^ ^ 

This gives two branches of dispersion curve. In homogeneous isotropic media, the two 
modes correspond to right and left circular polarizations. It should be noted that the 
condition for the existence of real solution is eju > 0. This means that the left-handed 
medium is defined in the same way as in non-chiral medium, in contrast to the case of 
ChC eqs. 


2.2 Case of ChC eqs 

A same way of solution is possible in this case, too. After eliminating £>, B from the 
ChC eqs, the solution for V x H, V x E has a same form as the one, where a, b, d, e of 
previous section are replaced with the following f,g,h,j, respectively 

/ = , (26) 
g = -i{uj/c){e + fiiri) , (27) 

h = —i{u/c)fi, (28) 

j = -i{uj/c)jlr] . (29) 

Further transformation of the equations of V x H, V x E into a set of homogeneous 
equations of E, H allows us to obtain the dispersion equation, as the condition for the 
existence of non-trivial solution, 

1 _ 

— = - 0 ± \/{A(^-0}^ + 4e/i ] , (30) 

oj z 

where we take all the combinations of ±. The condition for real solutions is 

+ > 0 , (31) 

which is less restrictive than non-chiral case. 


3 Discussions 

First of all, we note that, for both DBF and ChC eqs, the well-known result of [ckju'f' = 
e/r is obtained in the absence of chirality (/9 = 0 and .^ = 0, ?] = 0). Also both of the 
constitutive equations exhibit the typical behavior of chiral medium, i.e., the existence 
of the two branches with polarization dependent refractive indices. In the non-resonant 
region, both of them could be used to fit experimental results via appropriate choice of 
parameter values. 
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3.1 Resonant region of left-handed chiral medium 

A decisive difference appears in resonant region. An example will be the left-handed be¬ 
havior emerging in the neighborhood of a chiral resonance with El-Ml mixed character. 
Such a case has been treated by the first-principles theory of macroscopic constitutive 
equation iQ (Sec.3.8.1 and Sec.4.1.1). It shows a pair of dispersion curves for left and 
right circularly polarized modes, which have a linear crossing at fc = 0. It will be a 
test for the phenomenologies whether such a linear crossing can be realized or not by 
choosing parameter values. 

The dispersion equation in the first-principles macroscopic formalism is 

0 = det|^-l - [l-f —> (32) 

where Xe2(fc,a;) is the T component of susceptibility tensor (sec.2.5 of Q). Let us 

CT) 

choose a chiral form of susceptibility tensor Xem as 

1 + = (eb + a' + c'/c^)! + [ ] , (33) 

where the terms with a', 6', d represents the contribution of a pole ~ 1/{ojq — oj) with 
mixed (El, Ml) character, while Cb is the background dielectgric constant due to all 
the other resonances. We assume that the resonance with mixed El and Ml characters 
occurs in the frquency region of Cb < 0, i.e., a chiral version of left-handed medium. 
The dispersion equation is 

= ejl±j3jl — ^ (34) 

and its solution is given as 

_ j _ 

— = ±-[±/3/fi++ ^ (35) 

where we take all the combinations of ±, and 

(3 = ojb'/c , e = eh + a' ft = 1/[1 — {oj/cYd] . (36) 

Noting that h' is a chiral parameter corresponding to of the ChC eqs, we see that this 
equation is the same type as eq.(l30l). but not as eq.(l25l). An example of this dispersion 
relation is given in Eig.l. The characteristic behavior of the dispersion curves is a linear 
crossing at k=0 (Eig.4.1 of [171). 

In order to check whether the dispersion equations (l25l and (l30l) have this typical 
behavior of ’’linear crossing at /c = 0”, we focus on the behavior of the dispersion 
equations near fc = 0. Since both of the dispersion equations are given in the form 
"ckjoj = F{uy\ we only need to examine how the function F{u) approaches to zero in 
each case. 

The microscopic model of left-handed chiral medium given above consists of a mat¬ 
ter excitation level with (El, Ml) mixed character in the frequency range of Cb < 0. 
This means that all of XeE, XeB, XmE, XmB have a common pole (at u = Uq) and XeE is 
largely negative in the frequency range of interest to make Cb < 0. Namely, 
of ChC eqs have a common pole dX oj = ojq and e = Cb -f a /(ccq — oj). It may appear 
that the r.h.s. of eq. (l30l becoms zero for frequency satisfying /t = 0 or e = 0. However, 
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Figure 1: Dispersion curves of a chiral left-handed medium for the model in the text. 
Both ordinate and abscissa are normalized by the frequency of the pole as oj/ojq and 
ck/ojQ. Two horizontal lines show the frequncies of e = 0 and /r = 0. 


as mentioned before, the zero of /t corresponds to the pole of Ml transition, which in 
this case is common to the pole of ^ and t). Therefore the zero of fi is cancelled in the 
product — 7]). Thus, the only remaining possibility of zero arises from e = 0. The 
cu-dependence of the r.h.s. eq.dSO]) near zero point can be found by rewriting it as 

± i-. (37) 

2 A(? -'() ± ViMi-1)^+ 

At the frequency satisfying e = 0, i.e., + a' = 0, all of /i, r] remain finite, and one 

of the ± combinations in the denominator remains finite, so that the whole expression 
becomes zero for this combination. This occurs for both signs of ± in front of the whole 
expression. For negative Cb and positive numerator of the pole ~ 1 /(ccq — ce), Cb + a' = 0 
occurs at a; = cCz < ccq and in its neighborhood e ~ (ce — oj^). This shows that the r.h.s. 
of eq .(l30l) behaves like ~ (oj — Uz), which means the linear crossing of the two branches. 
Note also that uj^ lies inside the frequency range of eq.(l3T]). 

Now we check whether the same behavior is obtained for DBF eqs by assuming eq. 
dH), from which we obtain 


1 


1 

A 





(38) 


This shows that l//i has the same pole as l//i at cu = ojq, so that the factor /r//i on the 
r.h.s. of the second equation does not have the pole at ui = Uq via cancellation. There¬ 
fore, there is no chance for the denominator of the r.h.s. of eq.(l25]) to diverge. Hence 
the only possibility of its becoming zero comes from the factor y/ejl on the numerator. 
In view of the fact that the zeros of e, /r occur at different ui’s, e.g., at ui^i and U 22 (i^zi 
> U 22 ), the ce-dependence of should be ~ — oj or ~ y/ui — 10^2 in the neig- 

borhood of the zeros. Therefore, no linear crossing is possible in the DBF dispersion 
curves. The two zeros are the boundaries of the region of left-handed behavior. 
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One might argue that other type of ce-dependenee than eq. dH) eould lead to the linear 
erossing behavior. But one eannot freely give the ce-dependenee even as a phenomenol¬ 
ogy. Linear suseeptibilities should be a sum of single pole funetions. In the absenee of 
the first-prineiples theory for DBF eqs, it would be quite diffieult to give an appropriate 
model on a reliable basis. 


3.2 Conventionality vs. Logical Consistency 

DBF eqs have been popularly used in the maeroseopie argument of ehiral systems, es- 
peeially in the field of metamaterials researeh. As long as they are used for nonresonant 
phenomena as a praetieal tool, there is not mueh to say against it, exeept for the dif- 
fieulty in assigning mieroseopie meaning to the parameter jS. However, the restrietion 
to the nonresonant phenomena does not seem to be widely reeognized, to the knowl¬ 
edge of the present author. In faet there are examples of its use for resonant phenom¬ 
ena [fTTl [T 2 I. (The eonstitutive equations used in [fT^ are not exaetly DBF eqs, but 
D = eE — i^H , B = /xiT -|- i^E, different also from ChC eqs.) 

From the qualitative differenee of the two dispersion equations (l25l) and (l30l) in res¬ 
onant region, and from the faet that DBF eqs have no support by mieroseopie theory in 
contrast to ChC eqs, the use of DBF eqs for resonant phenomena is risky. As a conven¬ 
tional approach with a long history, DBF eqs might be kept in use further, but the validity 
limit should be kept in mind. However, if we consider that ChC eqs can be handled as 
easily as DBF eqs, and that they are consistent with the microscopically derived macro¬ 
scopic constitutive equationn, it is highly recommended to use ChC eqs. For problems 
requiring severe distinction, logical consistency should be preferred to conventionality. 


4 Conclusion 

The DBF eqs, popularly used as constitutive equations of chiral media, should be re¬ 
garded as a phenomenological theory applicable only in nonresonant region. In resonant 
region, it would lead to a qualitatively erroneous result. On the other hand, the ChC eqs, 
consistent with the first-principles microscopic constitutive equations, can be used both 
for resonant and nonresonant problems. 
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